Super-inflation in Loop Quantum Cosmology 
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We investigate the dynamics of super-inflation in two versions of Loop Quantum Cosmology, one 
in which the Friedmann equation is modified by the presence of inverse volume corrections, and 
one in which quadratic corrections are important. Computing the tilt of the power spectrum of the 
perturbed scalar field in terms of fast-roll parameters, we conclude that the first case leads to a 
power spectrum that is scale invariant for steep power law negative potentials and for the second 
case, scale invariance is obtained for positive potentials that asymptote to a constant value for large 
values of the scalar field. It is found that in both cases, the horizon problem is solved with only a 
few e-folds of super-inflationary evolution. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

An inflationary epoch is currently the most promising 
model for the origin of large-scale structure in the uni- 
verse [l| . The predictions of inflation are fully compatible 
with the most recent observations suggesting that struc- 
ture originated from a pattern of near scale-invariant, 
Gaussian and adiabatic primordial density fluctuations 
Despite these successes, however, a number of impor- 
tant questions remain. In particular, in what fundamen- 
tal theory will inflation be seen to arise? Having asked 
the question, it is worth noting that there have been suc- 
cessful implementations of inflation both in the context 
of ordinary field theory [3L H| as well as in the context of 
string and M-theory 0, |, |fl 0, 1, I, M, HO, El M Q ■ 
Given the importance of inflation and the need to ex- 
plore all possibilities of accommodating it in alternative 
theories of quantum gravity, in this paper we turn our 
attention to inflation in the context of Loop Quantum 
Gravity (LQG) [H. 

LQG is a background independent and non- 
perturbative canonical quantisation of general relativity 
based on Ashtekar variables: su(2) valued connections 
and conjugate triads. The variables used in the quanti- 
sation scheme are then holonomies of the connection and 
fluxes of the triad. The restriction of LQG to symmet- 
ric states gives rise to Loop Quantum Cosmology (LQC) 
p^ |. Although it is a particular limit of the more gen- 
eral LQG, and therefore can not be said to have generic 
features, LQC has produced a number of intriguing re- 
sults and resolved many problematic issues present in the 
earlier Wheeler de Witt quantum cosmology. In partic- 
ular, LQC can lead to a non-singular quantum evolution 
with the origin of the non-singular behaviour bc- 
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ing traced to the methods used to quantise inverse triad 
operators in LQG fTsLH^]. 

While the consequences of this quantum evolution are 
fascinating, it is difficult to connect it to existing theories 
of the early universe which tend to be based on classical 
dynamics, and in particular to inflation. Therefore, an 
approach to LQC has been developed in which effective 
or 'semi-classical' equations are derived and studied. In 
the isotropic setting, the effective equations which have 
been studied predominately to date include high energy 
modifications to the classical dynamics which originate 
from the spectra of q uantum operators related to the in- 
verse scale factor [I9l.l20l.l2ll.l22j]. In the context of scalar 
field driven inflation, a number of important effects follow 
from these modifications. These include the possibility 
that the field can be excited up its self-interaction po- 
tential [H, [H [H [13, leading to a subsequent 
period of slow-roll inflation. Most intriguing of all, how- 
ever, is the presence of a super-inflationary phase which 
occurs during the early phases of the universe's evolution 
irrespective of the form of the field's potential [111, [2§| . 

More recently, however, a further 'semi-classical' mod- 
ification, which arises from the use of holonomies as a ba- 
sic variable in the quantisation scheme, has been derived 
in the isotropic setting [301 ] . The modification is remark- 
ably simple and takes the form of an additional negative 
p term in the effective Friedmann equation, which ap- 
pears in addition to the usual positive p term. Such a 
term has a number of effects, it forces a collapsing uni- 
verse to undergo a non-singular bounce once a critical 
density is reached, and immediately after this bounce it 
also causes the universe to undergo a period of super- 
inflation. It is intriguing to note that such a term also 
appears in braneworld models with an extra time-like di- 
mension [3lj |. 

While the two sets of modifications discussed above 
have rather different origins, it appears that the quali- 
tative effects of both the inverse scale factor effects and 
the p 2 term are rather similar. In particular, they both 
give rise to a period of super-inflation during which the 
Hubble factor rapidly increases, rather than remaining 
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nearly constant as is the case during standard slow-roll 
inflation. Given that such a super-accelerating phase ap- 
pears to be a robust prediction of LQC, it is important 
to study both the background dynamics, and particularly 
the cosmological perturbations, which such a phase gives 
rise to. Considering a universe sourced by a scalar field, a 
number of important results have already been obtained. 
A scaling solution for the effective equations which arise 
from the inverse scale factor modifications has been de- 
rived [33 |. and a number of attempts made at stud yin g 
perturbations in the super- inflationary regime (33l . |34| . 
Similarly the scaling solution for the p 2 effective equa- 
tion has also been derived [H, Hi| . 

It is also interesting to note the close connections be- 
tween the super-inflationary phases in LQC and the evo- 
lution of a universe sourced by a phantom field, and 
with the e kpyr otic evolution of a collapsing universe 
[37L [H, [H, HO] . In all these cases the magnitude of 
the Hubble rate grows with time. Moreover, the scale 
factor duality discussed in [4lj maps the ekpyrotic col- 
lapse onto the super-inflationary scaling solution for the 
inverse scale factor modified equations. On the other 
hand, another duality maps the ekpyrotic collapse phase 
onto the dynamics of a universe sourced by a phantom 
field [4l|. These three regimes are therefore all related 
to one another. Furthermore, given that the ekpyrotic 
collapse is thought to offer a method for the generation 
of scale-invariant perturbations [42[ (as is the dual super- 
inflationary phase sourced by a phantom field [43] ) , it is 
reasonable to expect that a similar mechanism may oper- 
ate in the super-inflationary phases of LQC. Indeed such 
a mechanism has been discussed previously |34j |. though 
its relation to ekpyrotic and phantom models was not 
emphasised. 

In this study we aim to explore further the phe- 
nomenology of super-inflation in LQC. One complication, 
however, is that the relative status of the two sets of 
modifications discussed is at present unclear. We there- 
fore take a pragmatic approach and study the dynamics 
when each of the modifications is considered in turn, but 
not including both sets of modifications simultaneously, 
although we believe it should not be too difficult to in- 
corporate them both. A further difficulty is that despite 
considerable progress the understanding of metric pertur- 
bations in LQC is at present incomplete We there- 
fore restrict our attention to perturbations in the scalar 
field as a first approximation. This approach allows us to 
establish a framework for dealing with perturbations in 
LQC in which metric perturbations can be incorporated 
as our understanding advances. 

The paper is organised as follows. In section [TTJ we in- 
troduce the cosmological evolution equations which arise 
in LQC including the inverse volume corrections. Solu- 
tions are obtained including those showing scaling be- 
haviour, and the primordial spectrum of scalar perturba- 
tions is calculated for each solution in terms of 'fast-roll' 
parameters. The stability of these solutions is then dis- 
cussed. In section Inll we analyse the dynamics when the 



modification is induced by a p 2 correction to the Fried- 
mann equation. Concentrating on the evolution just after 
the bounce, we demonstrate the existence of the super- 
inflationary solution, obtain the scaling dynamics of the 
system, the primordial spectrum of scalar perturbations 
as well as the stability of the background solutions. Sec- 
tion llVl discusses the way in which super-inflation in LQC 
can solve the horizon problem with a small number of e- 
foldings and we conclude in section fVl 



II. EFFECTIVE FIELD EQUATIONS WITH LQC 
INVERSE VOLUME CORRECTIONS 



The first set of modified equations which we consider 
are those which incorporate two functions, Dj t i(a) and 
Sj(a), into the dynamics. These functions arise because 
of the presence of powers of the inverse scale factor in the 
Hamiltonian constraint for an isotropic and homogeneous 
universe. A full discussion of the origin of these terms can 
be found in Ref. [20;, 22] (a summary can be found in ap- 
pendix B of Ref. (ifl), but here we simply state their 
basic properties. We are implicitly considering cither 
positively curved or topologically compact models. This 
ensures that the size of the fiducial cell does not enter 
in the equations of motion. D and S are both functions 
of the scale factor, and their form changes depending on 
the values of two ambiguity parameters: I which takes 
values in the range < I < 1, and j which takes half 
integer values. When the scale factor approaches zero, D 
and S also approach zero, whereas as a increases above 
the critical value a* which depends on j, they both tend 
to unity. 

The modified Friedmann equation is given by 



H 



(1) 



where a dot denotes differentiation with respect to cosmic 
time t, and k 2 = 8irG. In what follows we choose units in 
which k = 1. We have omitted the curvature contribu- 
tion as we assume either a compact flat universe or that 
the curvature term rapidly becomes sub-dominant and 
that it can be safely neglected. The equation of motion 
for the scalar field takes the form 
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A subscript <j) means differentiation with respect to the 
field. These equations can also be combined to give the 
Raychaudhuri equation 
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A. Scaling dynamics 

We will be interested in the regime a <C a*, 
where the function Dji(a) may be approximated by 
a power law of the form D(a) = D+a n , with = 
(3/(3 + 2l)f /2{1 ~ l) o-Jd-sVCi-O and n = 3(3 _ _ /) 
takes values in the range 9 < n < 00. Likewise, the func- 
tion S(a) may be similarly approximated by S(a) = S l *a r , 
where S+ = (3/2)a~ 3 and r = 3, though we keep r ar- 
bitrary in our calculations for generality. For a 3> a* , 
S± ~ D* w 1 and r = n = 0. Inserting this form for the 
functions S and D into Eq. f3|, we can clearly see that 
for an expanding universe, and with n > 6 + r which oc- 
curs for all /, H is necessarily positive (assuming that the 
potential is either positive or the term involving SV can 
be neglected). Hence super- inflation is occurring. We 
will confine ourselves to these situations in what follows. 

To study this regime further, it proves convenient to 
introduce the variables 



(4) 



where p = cb 2 /2D + V(4>). Using these definitions, the 
equation of motion for the scalar field © can be written 
for an expanding universe in terms of a system of first 
order differential equations as 
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vv a 




V: 



(5) 

(6) 
(7) 

(8) 



and the third point is a scaling solution for which the 
kinetic and potential energies evolve in a constant ratio 
to one another: 



A 



r = i + ^Cn-r) 



(13) 



The scaling solution is therefore well defined for A 2 < 
6a 2 for positive potentials and for A 2 > 6a 2 for negative 
potentials. In the remainder of this analysis, we will focus 
mainly on the scaling solution for negative potentials as 
this is the case that, as we shall see, leads to a scale 
invariant power spectrum of the perturbed field. For this 
case one can check that the universe is undergoing super- 
inflationary expansion. 

Considering the fixed point for the scaling solution 
(fT2| . one can write 




A 



which upon integration gives 

2A 



\/6a 



(n — r)a 



(14) 



(15) 



where we have set the integration constant to zero with- 
out loss of generality. 

Then inserting this relation into the definition of A in 
© gives 



V = V 



A/3 



where f3 — — 2A 2 / (n — r)a > 0. 

Considering now the fixed point for y we have 



V 
P 
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3~ff2 



1 



_*1 

6a 2 



(16) 



(17) 



with a = 1 — n/6 < and N = In a. 
subject to the constraint equation 



2 1 2 
x ±y 



1. 



These variables are 



(9) 



The plus and minus signs correspond to positive and neg- 
ative potentials, respectively. Using the constraint equa- 
tion to substitute for y in Eq. §5§ renders Eq. © redun- 
dant. 

The resulting system defined by Eqs. ([5]) together with 
the constraint equation and ([7]) has three fixed points for 
A =/= 0. Two of them represent kinetic energy dominated 
solutions, valid for all values of A: 

x = -1, y = 0, r = l-^(n-r), (10) 
x = +1, y = 0, r = l + ^(n-r), (11) 



Differentiating Eq. (|15p and eliminating <f> using Eq. JTj 
then substituting for p in terms of V using Eq. Q17p and 
finally substituting for V in terms of a using Eq. (I16|) 
and Eq. (|15|) . one obtains an expression between a and 
d. Then by integrating this expression we can determine 
that the scale factor evolves as a power law in time while 
the universe evolves according to the scaling solution. In 
terms of conformal time dt — adr, we find 



a{r) = {-rf , 



(18) 



where for an expanding universe r is negative and in- 
creasing towards zero and 



2a 



P = 



21- (2 + r)a 



(19) 



where, for direct comparison with previous literature, we 
have introduced the slow-roll parameter e = A 2 /2 and 
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A = — v2e for <f> > 0. Using this form of a we find that 
H = p/dT, and it is straightforward to show that 



<f>'(r) 



V(r) 



2V2 



2e- (2 + r)a V 5 t 
4(3a 2 - e) 1 



(2e- (2 + r)a) 2 S(ar) 2 ' 



(20) 
(21) 



where a prime means differentiation with respect to con- 
formal time, r. Equations (fT8|) - (f2"Tj) form the basis of 
our analysis. This scaling relation (for S = 1) was first 
uncovered in Ref. [32] using a different procedure. 



B. Power spectrum of the perturbed field 

For a universe which evolves according to the scaling 
solution, the primordial spectrum of scalar perturbations 
produced by this super- inflationary phase was previously 
calculated in Ref. [34|. It was found that the spectrum 
tends to exact scale invariance for /3 > 1 (i.e. e > 1), 
without any fine tuning of the quantisation parameters of 
LQC. The purpose of this section is to review how scale 
invariance arises for the scaling solution with j3 3> 1, and 
to generalise the analysis of [34[ in order to allow for po- 
tentials which do not give rise to exact scaling solutions. 

In order to calculate the spectrum of perturbations, we 
now perturb the scalar field equation. The perturbation 
in the field 8(j> then satisfies the equation 



a 



D 



■D[V 2 -a 2 V^] 



which can be written in the form [341 ] 

u" + (-DV 2 + m 2 cS ) u = 0. 



(22) 



(23) 



where u is defined as u — aD 1 / 2 5(j) and the effective 
mass of the field u is given by 



(aD-^)" , 2 



'off 



DV, 



aD- 1 / 2 

Decomposing u in Fourier modes Wk, that satisfy 

wl+(-Dk 2 +ml s ) w k =0, 
the power spectrum is then given by 



(24) 



(25) 



(26) 



It was shown in Ref. [34( that the general solution to 
Eq. USD is 



Mr) = (4^2$(*) 

~H\l\{x)) (27) 



whenever m e gT is constant, where 

L ._ Vl-4mgffT3 
2 + np 



(28) 



and di and di are constants subject to the condition 
Mil 2 — M2I 2 = 1 and HPhx) and hW(x) are Hankel 
functions of the first and second kind, respectively. In 
the long wavelength limit, the power spectrum yields 



V u oc k 



3-2 \v\ 



l-\v\(np+2) 



(29) 



Scale invariance of the power spectrum is then attained 
when the spectral tilt An u = 3 — 2\v\ is zero. Since 
for a universe evolving according to the scaling solution 
Eqn. (TTHD 



■/?;., 



eV 2 



2+(3-2n)p+-(6 + 2n-n 2 )p 2 , (30) 



we can see from Eq. (|28p that scale invariance occurs 
whenever p — > 0, which, as we referred to, does indeed 
imply that e> 1 and consequently V < from Eq. (|21[) . 
There is one other value of p for which scale invariance is 
attained, p = — 4/(n + 4), however, we will not consider 
it any further. 

We would now like to generalise the form of the poten- 
tial we are dealing with so that it no longer has to be of 
exactly the form which gives rise to the scaling solution. 
For standard slow-roll inflation, where the kinetic energy 
is small compared with the potential energy, it is possi- 
ble to account for potentials of a form more general than 
a scaling potential by introducing slow-roll parameters. 
Such parameters parametrise the steepness of the poten- 
tial, and how this steepness evolves as the field moves 
along the potential. For a given field potential they also 
allow the dynamics which follow from a more general po- 
tential to be expanded locally about the dynamics which 
follow from a scaling potential with the same local steep- 
ness. The power spectrum which follows from the general 
potential can then also be written in terms of the slow- 
roll parameters. 

For the case at hand we would like to develop a similar 
expansion scheme. However, for the regime which we are 
considering in which e S> 1, it is clear that the kinetic 
energy is of approximately the same magnitude as the 
potential energy, and therefore the slow-roll approxima- 
tion is inadequate. Indeed in this case the field is evolv- 
ing rapidly along a steep negative potential, and we refer 
to the evolution as the 'fast-roll' regime. Our strategy 
will therefore be to determine other suitable small pa- 
rameters which characterise the steepness and curvature 
of the potential, and which we will refer to as 'fast-roll' 
parameters. The derived parameters we will arrive at 
are similar to those obtained in Ref. [12] , where fast-roll 
parameters were required to parametrise general poten- 
tials in the ekpyrotic scenario. This similarity is natural 
since, as we have already mentioned in the introduction, 
the evolution of the super-inflationary scaling solution in 
LQC is dual to the ekpyrotic collapse. 
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The first step in accommodating a more general class of 
potentials is to allow e to become time dependent. From 
its definition it then follows that 



-(2e) 3 / 2 W^', 



where we have defined 

V^V 1 V 



77=1 



2 V a 



D,4, 
D 



IT 



(31) 



(32) 



which can also be written in terms of the background 
quantities as 



, V M V 1 , n V a 1 



V 



(33) 



Likewise, we can calculate rj' in terms of a third param- 
eter £ 2 , 




(34) 



where 



VcfycfyVcf) 
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In particular, for the scaling solution where e is constant, 
one can verify that 77 = £ 2 = exactly. Since we are con- 
sidering potentials which are close to the form of a scaling 
potential, we expect that there will be solutions to the 
equations of motion of a form very similar to that given 
by Eqs. ([IS]). ([15]). (|2"0|) and ![2T]). when e is slowly vary- 
ing. Assuming that Eqs. (fl"8|) , (fT9j) and ([20f are indeed 
good approximations, we have in general that 



dine 
din a 



4 £_, ^«2i£. (36) 

a a In a 77 a 



Imposing that e and 77 are slowly varying, and since e is 
large in the regime which gives rise to scale invariance, 
requires 77 and £ 2 to be small, i.e. the potential must 
be nearly power law in form which is in agreement with 
our assumption. We refer to 77 and £ as the second and 
third fast-roll parameters, and for convenience introduce 
the first fast-roll parameter as e = l/2e, where in terms 
of relative magnitude we find e ~ 77 and £ 2 ~ C(e 2 ). 
Using these relations it is possible to verify by substitut- 
ing Eqs. (|18p - (|2"Tj) into the Friedmann and equation of 
motion that these solutions are indeed valid up to sec- 
ond order in fast-roll parameters for a general negative 
potential confirming the consistency of our analysis. 



Then using Eqs. (p~8|) - (f2T|) to substitute for the respec- 
tive quantities in the effective mass (I24[) , and expanding 
to first order in the parameters e and 77, we obtain 



4c 



n ( n \ r 

1 (H r) - - 

12 V 6 J 2. 



477 . 



(37) 



We note that although we have used the solution to 
Eq. (|2l)|) which is valid only when m G ffr is constant, the 
solution will remain sufficiently accurate provided that 
e does not vary significantly as a given k mode crosses 
the horizon from the small wavelength to the long wave- 
length regime. This is simply the condition that 77 is 
small, which we have assumed already. The spectral tilt 
can therefore be calculated at any given scale by insert- 
ing into Eq. ([3T)l the values that e and 77 take as this 
scale crosses the horizon. In particular, to compare with 
observations we need to consider the mode which cor- 
responds to the largest scales on the Cosmic Microwave 
Background (CMB). 

Finally we also note that we could have included time 
derivatives of the power n and r. The generalisation to do 
so is straightforward but for simplicity the computation 
we have presented does not include this possibility. 



C. Stability of the fixed points 

The analysis we have performed so far is intriguing. It 
appears that in the a <C a* regime, the scaling solution 
which follows from a steep negative potential can give rise 
to a scale-invariant power spectrum of scalar field pertur- 
bations, and moreover we can generalise the analysis to 
potentials which deviate from the scaling potential. 

However, there is another element which is involved 
in building a convincing theory for the origin of scale- 
invariant perturbations. That is, it would be highly de- 
sirable if the scaling solution was an attractor, so that 
initial conditions not exactly on the solution would evolve 
towards it, and the solution's stability against small local 
perturbations would be assured. 

To determine the stability of the scaling solution, we 
study the nature of the fixed points of Eqs. ([5]) and 
([7]). We do this by linearising the equations about the 
fixed points and determining the corresponding eigenval- 
ues (uj) in each case. For the kinetic energy dominated 
solutions, valid for an arbitrary A, in Eqs. (fT0|) and (jlip . 
we find their respective eigenvalues to be 



lu + = 6a — Vh~\ , 



w _=-i(n-r), (38) 
w _ = -i(n-r). (39) 



Since n > r and a < 0, the first fixed point is stable for 
A < — \/6a and the second for A > ^/6a. Turning to the 
scaling solution, Eq. (p~2|) . we find 



(40) 



u± = -— (9 ± J9 2 + 8a{n-r)(X 2 -6a 2 ) 
4a V / 
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V = V o 2 

81 — 1 
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1 0.5 1 1.5 

lna/a* 

FIG. 1: The evolution of the ratio —(j) 2 /2DV obtained by nu- 
merically solving the equations of motion for three different 
initial conditions (solid line). They approach the scaling so- 
lution given by x 2 /y 2 where x and y are given by Eqs. fjl2p 
(dashed line). We used as parameters Vo = — 10 -20 , 0i n it = 1, 
in Planck units and n — 15, r — 3 and aimt = 0.9a*, in a flat 
universe. 



Once again we are assuming either a flat universe or that 
the curvature contribution can be safely neglected. It is 
interesting that this form of the Friedmann equation is 
identical to the form which arises in braneworld scenarios 
with a single time-like extra dimension in the absence of a 
black hole in the bulk spacetime [31] . In the braneworld 
case a represents the brane tension while for the LQC 
case 2a represents the critical energy density arising from 
quantum geometry effects which leads to the scale factor 
undergoing a bounce as p approaches it. 

We are interested in high density regimes where p ap- 
proaches the bounding value of 2a. In this case, the 
term within brackets tends to zero, and the behaviour 
of the equations alters significantly compared with the 
classical behaviour. Indeed in this regime we have H > 
and for an expanding universe super-inflation takes place. 
We will again consider a scalar field dominated universe, 
hence, p — (f) 2 /2 + V((f>). We stress that we are study- 
ing inverse volume and quadratic corrections separately, 
hence we do not include the D and S functions in the 
Friedmann equation and in the definition of energy den- 
sity. The scalar field equation of motion 

4> + 3tf </> + V.4, = , (42) 

is unchanged from the classical form. 



where 6 = a (6 — r) — A 2 < 0. The scaling solution 
is therefore stable whenever A 2 > 6a 2 which coincides 
with the region of existence of this solution. In Fig. [T] 
we show the evolution of the ratio ~cf> 2 /2DV obtained 
by numerically solving the equations of motion. We can 
see that it approaches the value given by x 2 /y 2 where x 
and y are given by Eqs. (|12[) . Tyically the evolution only 
reaches the attractor when a > a* which is far outside 
the region where the appoximation D « a n is valid. We 
conclude that this solution must be extremely fine tuned 
in order to deliver the dynamics and power spectrum as 
described in the previous subsections. 

In the second part of this article we will be dealing 
with a second possibility of obtaining a super-inflationary 
regime that also leads to a scale invariant power spec- 
trum with the advantage that the stability of the scaling 
solution is no longer a dangerous issue. 



A. Scaling dynamics 

It was shown in Ref. [36| (also see [Hj]) that the effec- 
tive equations ([4"l"j) - ([4*2"j) also allow a scaling solution in 
which the kinetic and potential energy vary in propor- 
tion to one another. In this case the potential must be 
of the form V = Vbcosh(</>). However, while earlier works 
claimed otherwise, it can be shown that the scaling solu- 
tion is not an attractor. Moreover, the scaling solution 
which was found implies an evolution for the scale factor 
which is not of a power law kind, and this means that 
it is unlikely to give rise to a scale-invariant spectrum of 
perturbations. Given these difficulties, instead of focus- 
ing on the scaling solution we will ask whether there is a 
form of the scalar field potential which does result in a 
power law evolution of the scale factor. We are interested 
in the regime in which p ss 2a, when super-inflation oc- 
curs, and where H w 0. Inserting the power-law ansatz, 



III. EFFECTIVE DYNAMICS WITH 
QUADRATIC CORRECTIONS 

The second modification to classical dynamics which 
we consider follows from considering that holonomies are 
the basic variables for quantisation in LQC. This modifi- 
cation gives rise to a Friedmann equation of the following 
form 13011: 



H 2 



1 



2a 



(41) 



o(t) = H) ro , (43) 
where m < 0, into the time derivative of the Hubble rate, 

(44) 



2 V a 



we see that for p w 2a the kinetic energy of the field is 
4> 2 /2 m —m/t 2 which upon integration gives 



4> w ±V = 2m Int. 



(45) 
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Using Eq. (|4Tj) and expanding in 6H 2 /a < 1 we have 
that 



2a - 3H Z 



(46) 



and it follows from the definition of the energy density of 
the field and Eq. ((45]) that 



V w 2a - U e 



-\<t> 



(47) 



where Uq = 3m 2 — m > and A = 1 / V— 2m. It is evident 
that, in this regime, scaling exists between V — 2a and 
the kinetic energy <j) 2 /2. We now look for a more precise 
description of the dynamics. The form of the potential 
(l47l) motivates us to define the new variables: 



u 



V2a~~^ 



(48) 



such that p < 2<T and V(<^>) = 2(7 — U(<fi). In terms of a 
system of first order differential equations, the equation 
of motion of the scalar field now reads, 



y,N 

A ,jv 



-3a; 



3x d 



- Aa;y - 3a; y 



-V6A 2 (r- l)x + 3a; 2 1 

P 



2a 



where A and T are defined as 
P 



A : 



' U 



uu 



<t><t> 



u z 



(49) 
(50) 

7 ' (51) 

(52) 



The variables x, and y are also related by the constraint 
condition 



2 2 i 

x -y = -1. 



(53) 



Considering the regime discussed above where 2a/ p rs 
1, we can see A is a constant and by integrating A, we see 
that the U part of the scalar potential is given by 



(54) 



as we were expecting from Eq. (j47]) . 

We consider the section of the phase space in which 
x < 0, y > 0, and A > 0. Taking A to be constant, and 
substituting the constraint equation (|53|) into Eq. (|49]) . 
results in an autonomous system with three fixed points. 
Two of them are non-physical solutions with 



x = ±i, y = 0, 
and the third is a scaling solution with 



A 



y 



i 



A 2 
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(55) 



(56) 



The scaling solution is valid for all real values of A. 

As in the previous case, it is straightforward to show 
that as the universe evolves according to this solution, 
the scale factor undergoes a power law evolution 



a(T) = (-r)*\ 



where 



(57) 



(58) 



and e is here defined as e = (U^/U) 2 /2 w A 2 /2. This 
is of course what we expect since we began by searching 
for such a solution using the ansatz Eq. (|43|) . The time 
derivative of the field and the potential yields, 



V2e 1 



V = 2a 



e + 1 r ' 

3 + 1 



(1 + e) 2 (ar) 2 ' 



(59) 
(60) 



We are now ready to compute the spectrum of the scalar 
field perturbations produced by this power-law solution. 



B. Power spectrum of the perturbed held 

In this section we follow the same approach we took in 
the previous analysis of the scalar field perturbations. In 
terms of conformal time, the perturbation equation for 
the scalar field, </>, can be written as 



" = -2-5<P' + (V 2 - a 2 V M ) • 



(61) 



which in turn can be written in terms of u = aS<p as 

u" + (-V 2 + m 2 ff ) u = Q , 

and the effective mass of the field u is 

a" 



(62) 



™cff 



+ a 2 V 



<t><t> 



(63) 



Decomposing u in Fourier modes Wk, that satisfy 

< + (- fc2 +™off) w k = Q, 
the power spectrum is then given by 



(64) 



(65) 



The general solution to Eq. ([23]) is 
v>*(t) = (d^H^ix) +d 2 ^H\l}{x)) ,(66) 
where the subscript \v\ is 

(67) 



^l-4m 2 cS T 2 
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and d\ and c?2 are constants subject to the condition 
1^1 2 _ |d 2 |2 _ i and H\l\{x) and h\1\{x) are Hankel 
functions of the first and second kind, respectively. For 
large wavelength modes, the power spectrum can be ap- 
proximated by 



l—2\v\ 



(68) 



By substituting Eqs. IJ57J - §U§ into Eq. ||63J), we find 



-2 + 3p + 3p 2 



(69) 



Comparing this case to our previous results, we expect to 
have scale invariance for p — > 0. Once again, therefore, 
scale invariance occurs when the field <fi is rolling down 
a steep potential and the kinetic energy is not negligible, 
but comparable to V — 2a. Hence the evolution should 
again be understood as a fast-roll regime. 

Clearly, scale invariance is also obtained for p — > — 1 or 
e <C 1 which corresponds to the standard slow-roll regime 
that we are not concerned with for the purposes of this 
work. 

In order to extend this analysis to general potentials, 
as we did for the previous system of modified equations 
we studied, we allow e to depend on conformal time such 
that 



e -' = -( 2e -) 3 /W, 

which defines the fast-roll parameter 77, 

U M U 



77=1- 



U 2 



(70) 



(71) 



Similarly, 77' can be written in terms of the next order 
fast-roll parameter £ 2 as 



77 



2 aJ 



(72) 



where we have again defined a further fast-roll parameter 
by e = l/2e. Again the assumption that m^T is nearly 
constant as a given mode evolves outside the cosmological 
horizon is valid. Hence the spectral tilt for a given k 
mode can be calculated by inserting the values the fast- 
roll parameters took as the mode crossed the horizon in 
Eq. (|75p . It is clear that the system we have investigated 
here results in a scale-invariant spectrum of scalar field 
perturbations for e< 1 and 77 <C 1. 



C. Stability of the fixed points 

Linearising the system (|4"9")) . using Eq. (T53")) about the 
fixed points we find the following eigenvalues 10: for the 
unphysical kinetic energy dominated solution, Eq. (|55[) . 



= 6 ± iV6A . 



(76) 



and hence this solution is unstable. While for the scaling 
solution, Eq. (|56|) 



1 



(6 + A 2 ), 



(77) 



with 



and hence this point is a stable attractor for all values of 
A. A numerical analysis of this system shown in Fig. [2] 
supports our analytical results presented here. The fig- 
ure shows the evolution of the ratio cb 2 /(2a — V) ob- 
tained by numerically solving the equations of motion 
for three different initial conditions. They approach the 
scaling solution given by 2a; 2 /y 2 where x and y are given 
by Eqs. ([56")) . When a ^> a* the quadratic corrections 
become negligeble and the numerical evolution diverges 
from this attractor. 



IV. NUMBER OF e-FOLDS 



e 2 = 



2 u^L)u^ (73) 



u 



u 2 



For the scaling solution, it can be verified that both 
77 and £ 2 vanish. As in the previous case, we can use 
Eqs. (|57|) , (|58|) and ([59)) as approximate solutions, and 
hence we have 



dine 
din a 



4e7?, 



<iln a 



n 



(74) 



which means that for a large and slowly varying e the 
parameter 77 must be small, and for a slowly varying 77 
the parameter £ 2 must also be small. Hence, the U part 
of the scalar potential must be close to exponential. 

Using now Eqs. (|F7|) - ([59")) in the expression for the ef- 
fective mass of u Eq. (f6"3"]) . where e is now time dependent, 
and then using Eq. (f6"T)) . we find that An„ = 3 — 2\u\ is 
to first order 



An,, 



-4(6-77) 



(75) 



Before concluding, it is important to address the ques- 
tion of whether a sufficient amount of super-inflation can 
occur in LQC to account for the largest scale perturba- 
tions observed on the CMB. This in turn is equivalent to 
asking whether the super-inflationary phases can solve 
the cosmological horizon problem. 

It is clear that only a small number of e-folds of super- 
inflation can be considered generic for either of the modi- 
fied sets of evolution equations we have studied [29j] . This 
might be considered disappointing, since experience from 
standard inflation suggests that approximately 60 e-folds 
of inflation are required for consistency with observa- 
tions. However, the inflationary periods we have been 
studying are considerably different to standard inflation 
and this has a dramatic effect. 

Solving the horizon problem is essentially the require- 
ment that aH grows sufficiently during an early stage 
of the universe's evolution. While in standard inflation 
this is accomplished by a changing rapidly as H remains 
nearly constant, in our case the converse appears to be 
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U = Uq exp(- 





hxa/ainit 

FIG. 2: The evolution of the ratio 4> 2 /(2a - V) obtained by 
numerically solving the equations of motion for three differ- 
ent initial conditions (solid lines). They approach the scaling 
solution given by 2x 2 /y 2 where x and y are given by Eqs. (|56p 
(dashed line). We used as parameters Uo = 10 -2 , cpimt = 1 
and a — 0.41 in Planck units, in a flat universe. 



possible, that H increases sufficiently as a remains nearly- 
constant. The number of e-folds usually only refers to 
the change in a, and so we only expect a small number 
of e-folds to be necessary in a super-inflationary phase, 
provided that H changes sufficiently. To confirm the ex- 
pectation that our super-inflationary phases can indeed 
solve the horizon problem, let us now quantify our qual- 
itative arguments. 

During super-inflation, perturbation modes exit the 
cosmological horizon, and once super-inflation ends, 
modes start to re-enter. Let us consider a perturba- 
tion mode with wavenumber k, such that k exited the 
cosmological horizon N(k) e-folds before the end of the 
super-inflationary phase, and re-entered sometime later. 
The mode re-entering the horizon today, fc* , must satisfy 
fe* = agHo, where subscript indicates quantities at the 
present epoch. Comparing this with the generic k mode 
we have: 



akH k 



QkHk a en d a re h Cleq -ffend 



ao-ffo ao-ffo Oend-Hend a rc h acq ao Hq 



(78) 



where subscript 'end' labels quantities at the end of in- 
flation, 'reh' at reheating, and 'eq' at matter radiation 
equality. Then employing the known evolution of the uni- 
verse from reheating to the present day, together with the 
measured value of the Hubble rate at the present epoch, 
Hq, and for simplicity assuming that the universe be- 
haves as if it is matter dominated between the end of 
inflation and reheating, we find: 



In 



ciqHo 



68 + In 



akHk 

^end-^end 



iln 
2 



Mi 



pi 



H, 



end 



in 



Pond 
. Preh . 



1/4 



(79) 



The energy scale at the end of inflation must be deter- 
mined by requiring that the magnitude of the curvature 
perturbation accounts for the temperature anisotropies in 
the CMB. Since we have only worked with the scalar field 
perturbation this is so far undetermined in our model and 
we therefore take -ff cn d to be the highest possible scale, 
i.e. H cn d = Mp\. A lower scale would lead to fewer re- 
quired e-folds. Further considering k = k* and assuming 
instant reheating, we find that 



hi 



(^end-"end 
ak, H k , 



68. 



(80) 



In order to determine the number of e-folds of super- 
inflation required, we now consider the cases in which 
the scale factor is undergoing pure power-law behaviour, 
a cx (— t) p . Using this together with Eq. ([50]) . we find 



and in turn 



JV(M = In 



^end 

a-k, 



68, 



= -68p. 



(81) 



(82) 



Recalling that p must be small and negative for scale 
invariance, we see that only a small number of e-folds are 
required. Although considering behaviour which deviates 
from pure power law behaviour will alter Eq. (|82p . it is 
clear that the conclusion of only a small number of e-folds 
being necessary will remain valid. 



V. DISCUSSION 

In this paper we have investigated the nature of super- 
inflation in Loop Quantum Cosmology. Considering 
two specific examples which lead to modifications of the 
Friedmann equation, one where the modification is due 
to the presence of inverse volume corrections, and the 
second where it is induced by the use of holonomies as 
the basic variable in the quantisation of Loop Quan- 
tum Gravity, we have demonstrated explicitly in both 
cases the existence of super-inflation solutions defined by 
H > where H is the Hubble parameter. Through the 
use of phase plane analysis we have been able to dis- 
cuss the nature of the attractor solutions and their sta- 
bility in both cases. Further we determined the scalar 
perturbations arising in these models, and through the 
introduction of fast roll parameters, showed that a class 
of potentials exhibit near scale invariant spectra. To be 
more specific we have shown that for the case with in- 
verse volume corrections, if we concentrate on the regime 
D a a™, three solutions exist, two of them correspond- 
ing to kinetic energy dominated solutions and one to a 
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super-inflationary scaling solution where the kinetic en- 
ergy scales in proportion to the potential energy. In this 
last case, when the potential is negative, scaling occurs 
when it is polynomial in form, and the perturbed scalar 
field equations for this potential are scale invariant (as 
first demonstrated in [34j). However, we have been able 
to go further. By introducing fast-roll parameters e and 
r\ in an analogous manner to the way introduced for the 
ekpyrotic scenario [42| . we were able to extend our cal- 
culation of the scalar power spectrum beyond the case of 
the polynomial scaling potential. In other words we can 
determine the degree to which scale invariance is broken 
in terms of the fast roll parameters. Through the stability 
analysis of these solutions we have seen that the scaling 
solution which gives rise to the scale invariant spectrum 
is a stable attractor. In general, however, the attractor 
is reached when a 3> a* when D « 1 and therefore we 
concluded that in only a limited range of parameters we 
have an analytical understanding of the dynamics of the 
system in the semi-classical phase a ^ a+. 

We considered a second set of corrections to the Fried- 
mann equation, arising from the quantisation procedure 
in LQG where a p 2 term in the Friedmann equation, anal- 
ogous to that found in particular braneworld models 31 1 , 
is present. In this case the scaling solution found in [361 ] . 
which involves a cosh(</>) potential, is not an attractor (in 
contradiction to the claims in Refs. [H, [36|) and we be- 
lieve it is unlikely to lead to a scale invariant spectrum. In 
order to obtain the solution with a scale invariant spec- 
trum we decided to take another route, namely search 
for the potential with the correct time evolution in the 
scale factor which would lead to scale invariance. In these 
models super-inflation occurs just after the bounce and 
using this fact we quickly arrive at the intriguing result 
that the corresponding potential is an uplifted negative 
exponential potential, and in this case the scaling solu- 
tion is stable. Scale invariance easily follows as does the 
generalisation to include a new set of fast roll param- 
eters, hence obtaining potentials which will yield small 
deviations from scale invariance. 

In both cases, we note that we do not expect the scalar 
potential to remain of the form we have used to give scale 
invariance throughout the entire evolution of the universe 
as a negative value of the potential after the end of super- 
inflation may lead to a subsequent recollapse. The full 
details of the exit from super-inflation, the subsequent 
form of the potential, as well as the transition to the ra- 
diation epoch, are beyond the scope of the present work. 

We emphasize that we only dealt with the evolution of 
the universe in the expanding phase, however, for closed 
models, the universe may have gone through a bounce 
from a collapsing phase into the super-inflationary evo- 
lution. It is natural to question whether the bounce is 



symmetric or not, i.e. if the period of super-inflation is 
nothing more than the counterpart of an identical defla- 
tionary period and an epoch of standard inflation is still 
necessary. Though this is a possible situation, we stress 
that the scaling solutions produced by the potentials we 
considered here are the attractors only during the ex- 
panding epoch, hence, the evolution through the bounce 
is indeed expected to be asymmetric. 

Finally, we have discussed the horizon problem in this 
set up. We found that for the class of potentials studied 
that lead to scale invariance, the required number of e- 
folds of super-inflation is of only a few, more specifically, 
N w — 68p where — 1 <C p < defines the time depen- 
dence of the scale factor of the universe, a = {—t) p . A 
previous criticism of infation within the inverse volume 
corrections approach has been that a* must be unphysi- 
cally large to give rise to 60 e-folds [46[ . We have seen in 
this work, however, that only a small number of e-folds 
of super-inflation are required, which suggests that this 
bound on a* can be evaded. We also note that previ- 
ous studies found a small probability of sufficient stan- 
dard inflation in the context of LQC [47]. However, their 
conclusions do not apply directly to our model of super- 
inflation since, as we have discussed, only a small number 
of e-folds are required and moreover, in their study, the 
Hubble rate was assumed to be constant. 

In our calculation we have neglected metric pertur- 
bations, instead we have only considered scalar pertur- 
bations and been able to find scale invariance through 
the scalar field perturbations. Given the similarity men- 
tioned throughout the paper between our LQC results 
and the ekpyrotic mechanism, it seems likely that when 
metric perturbations are included in the LQC calcula- 
tion, we will find that the scale invariant spectrum occurs 
in the Newtonian potential and the scalar field perturba- 
tion, but not in the curvature perturbation, just like the 
ekpyrotic case (for example see [42|), If this turns out to 
be the case, it would be interesting to consider whether 
the inclusion of a second scalar field in LQC, as in the 
'new ekpyrotic' mechanism, would allow the scale invari- 
ant spectrum to be transferred from the scalar field per- 
turbation into the curvature perturbation via an entropy 
perturbation |H, E • 
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